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ABSTRACT 

Partial frequency redistribution (PRD), describing the formation of the hne profile, has negligible 
observational effects for optical depths smaller than ~10^, at the resolving power of most current 
instruments. However, when the spectral resolution is sufficiently high, PRD modeling becomes es- 
sential in interpreting the line shapes and determining the total line fluxes. We demonstrate the 
effects of PRD on the H2 line profiles observed at high spectral resolution by the Far — Ultraviolet 
Spectroscopic Explorer (FUSE) in the atmospheres of Jupiter and Saturn. In these spectra, the 
asymmetric shapes of the lines in the Lyman (v"- 6) progression pumped by the solar Ly/3 are ex- 
plained by coherent scattering of the photons in the line wings. We introduce a simple computational 
approximation to mitigate the numerical difficulties of radiative transfer with PRD, and show that 
it reproduces the exact radiative transfer solution to better than 10%. The lines predicted by our 
radiative transfer model with PRD, including the H2 density and temperature distribution as a func- 
tion of height in the atmosphere, are in agreement with the line profiles observed by FU SE. We 
discuss the observational consequences of PRD, and show that this computational method also allows 
us to include PRD in modeling the continuum pumped H2 fluorescence, treating about 4000 lines 
simultaneously. 

Subject headings: planets and satellites: individual (Jupiter, Saturn) — ISM: molecules — ultraviolet: 
solar system - ultraviolet: ISM - radiative transfer — line: profiles 
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1. INTRODUCTION 

The study of the ultraviolet spectra of giant planets started with the Voyager fly-bys and rocket observations 
(jBroadfoot et al.lll977t iMoos et al.lll969| ). and continued with Galileo and Cassini, as well as with observations from 
space-based observatories such as the International Ultraviolet Explorer (lUE), the Hopkins Ultraviolet Telescope 
(HUT), and the Far-Ultraviolet Spectroscopic Explorer (FUSE). The increasing spectral resolution of the cited 
observations allowed for a better separation of the H2 excitation processes, constraining the thermal structure of the 
planetary atmospheres. The far-UV spec tra of giant planet atmospheres are dominated by molecular hydrogen emission 
produced by electron- impact excitation (|Liu fc Dalgarnolll996l : IWolven fc Feldnianlll998[ ). and, to a lesser extent, by 
fluorescent pumping (jFeldman et alJll993h . The molecular hydrogen transitions observed in the far-UV occur via 
de-excitation from the excited electronic states, with most of the contributions coming from the states B^E+. C^H„, 
EF^E+, B'^E+, D^H„, B"^E+, and D'^H„. Since these electronic states lie at 11 eV or more above the ground electronic 
state (X^E^), only photoelectrons, auroral electrons, and direct solar UV radiation have enough energy to excite these 
transitions. The solar far-UV continuum is however too weak to make a significant contribution to the observed 
spectrum, and the only observed radiatively-excited features arise from selective pumping of H2 transitions coincident 
with str ong solar emission lines . At H2 rotational te mperatures below 1000 K that characterize the atmospheres of giant 
planets ()Liu fc Dalgarnolll996t fWolven et al.lll997D . the main mechanism generating H2 fluorescent lines is pumping 
by solar Lyman-/?, which generates lines in the Lyman system (B^Ej-X^E^). Several of these lines were observed in 
the HUT observations of the Jovian dayglow by iFeldman et al.l (jl993D . In addi tion, weak fluorescent emission in these 
lines has been use d to detect H2 for th e first time in the atmosphere of Mars (jKrasnopolskv fc Feldmanll200lD and in 
a cometary coma (jFeldman et al.|[2002D . 

We report FUSE observations of Jupiter and Saturn, in the spectral range 905-1187 A at a spectral resolution 
of ~0. 08-0. 12 A, sufficient to resolve the line profiles of the H2 fluorescent lines pumped by solar Ly/3. For both 
planets, we observe a broad asymmetric profile for three of the Lyman (6— v")P(l) lines, namely Ly(6— 3)P(1) A1166.8, 
Ly(6— l)P(l) A1071.6, and Ly(6— 2)P(1) A1118.6, not reproduced by frequency-independent fluorescence models. These 
line profiles can be explained by the asymmetry in the overlap between the Ly(6— O)P(l) line and the solar Ly/3, 
combined with the coherent scattering of photons in the line wings. We treat separately the electron-impact and 
radiative excitation contribution to the spectra, and investigate both the complete frequency redistribution (CRD) 
approximation, as well as the PRD solution for the formation of the Ly/3-pumped H2 lines. For the line-pumping 
model, we build synthetic spectra starting from depth dependent atmospheric models, and show that the observed 
line profiles cannot be reproduced under the approximation of CRD, even taking into account the spatial variation of 
the atmospheric composition and temperature, or the differences in the numerical integration method for the radiative 
transfer equations. The use of PRD on the other hand correctly describes the coherent scattering in the line wings, 
proven by the agreement between the model and the data, and allows an independent confirmation of the H2 abundance 
and atmospheric structure. 

The P RD effects have been thoroughly studied for the formatio n of resonant and subordinate lines in stellar atmo- 
spheres (jShine et al.lll975l : iHubenv fc Liteslll995l : IUitenbroekll200l| ) , but only thcprofiles of resonant lines with overlap- 
ping transitions have been taken into c onsideration for planetary atmospheres ([Gladstone 1982 : Wallace fc Yellc, 19891 : 
lGrifiioenll200(3l : iBarthelemv et al.ll2004[ ). The importance of PRD has also been pointed out for the formation atomi c 
Hnes in quasi-stellar objects and active galactic n uclei (AGN s) (|Collin-Souffrin fc Dumontlll986l : [Avrett fc Loeseill988[ ). 
and for the X-ray Bowen fluorescence in AGNs (|Sakcill2003l ). The formalism presented in this paper is geared towards 
the treatment of cross-redistribution in fluorescent molecular lines, which has not been investigated to date. The case 
of molecular lines is complicated due to the multitude of transitions, the large wavelength range spanned by all possible 
subordinate lines, and the fine frequency grid required to accurately describe the line profiles. 

The FU SE observations of the atmospheres of Jupiter and Saturn represent the first detection of PRD effects in 
subordinate molecular fluorescent lines, and offer grounds for investigating the same process in other astrophysical 
environments. The far-UV H2 emission can be difficult to interpret due to the contributions of different excitation 
mechanisms, dust extinction and scattering, as well as the large optical depth effects at high column de nsities. Complex 
models of m olecular hy drogen emission i n astrop hysical environ ments have been constructed (e.g., [Sternberg 119891 : 
iDraine fc Bertoldi 199(i iLiu fc Dalgarna|[l996. : ,Shaw et al1l2005f ). taking into account the formation and destruction 
processes, the chemistry, and the detailed balance equations. However, in these models the radiative transfer of the 
molecular lines is often treated in a simplified manner, due to computational limitations, and does not attempt to 
reproduce the line profiles. Although the use of the PRD functions to characterize the emission line profiles involves 
adding an additional dimension to the radiative transfer problem and requires large computational resources, the 
consequences for the predicted spectrum can be significant. 

Our approach for modeling the H2 fluorescence is complementary to previous models, focusing on reproducing the 
line profiles, while making simplifying assumptions for the ground state population and the geometry of the system. 
In order to mitigate the computational difficulties, we investigate approximate solutions for the r adiative transfer 
equat ion, in the case of a plane-parallel atmosphere. Starting from an optically thin approximation (iLiu fc Dalgarnol 
|1996( ). we find that treating a medium of large optical depth as a series of optically thin layers reproduces the Feautrier 
solution for a single line at the 10% level. This multi-layer approach is the method of choice to generate the synthetic 
spectra for Ly/3 pumped fluorescence of H2, reproducing the FUSE observations. This method also allows us to test 
the effects of PRD on continuum pumped H2 fluorescence, a case involving several thousand transitions, that would 
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otherwise be intractable. In a simple toy model for the continuum pumped H2 fluorescence, using just a single layer 
approximation, the inclusion of PRD results in significant blending of the overlapping line wings, and has the effect of 
reducing the line contrast by generating a pseudo-continuum. We show that this process has important consequences 
for the detectability of lines for large optical depths and high spectral resolution observations. 

The FU SE observations are presented in Section [21 Section [3] describes the H2 fluorescence model, the approx- 
imations used, and the introduction of PRD to reproduce the line profiles. The application to line-pumped H2 in 
planetary atmospheres and the continuum pumping test case follow in Section 21 A summary of our conclusions is 
given in Section [5l 

2. FUSE OBSERVATIONS OF JUPITER AND SATURN 
2.1. Observations and Data Reduction 

The Jovian dayglow was recorded during drift scans of the eastern limb of Jupiter made on 2001 October 19 and 21 
using the FUSE MDRS (4" x 20") aperture. This approach was chosen to minimize the possibility of contamination 
by the much stronger polar aurora and to take advantage of possible limb brightening of D I Lyman-/3 emission. A 
single scan of ~ 3000 s was made during each of five consecutive orbits on both days. For each scan, the long side 
of the slit was aligned nearly north-south, along a chord centered on the Jovian equator. The scans started with the 
center of the slit positioned at 0.77 Rjup from the center of Jupiter, and continued towards the eastern limb a rate of 
~ 0'.'23 minute"^. The leading edge of the slit crossed the limb of the planet ~ 700 s after the start of the scan, and 
~ 300 s later the spacecraft crossed the terminator into orbit day. A montage of slit positions during these observations 
is shown in Figure [H At the time of the observations, Jupiter's heliocentric distance was 5.14 AU, and its geocentric 
distance was 4.86 and 4.83 AU for the two days, respectively. 

FUSE consists of four separate telescope/spectrograph assemblies, two employing lithium fluoride coatings (LiF) 
and two using silicon carbide (SiC). Each of four separate detectors (denoted la, lb, 2a, and 2b) records two spectra, 
one each from a LiF and a SiC channel spanning a g iven wavelength in terval , giving a tota l of eig ht separate spectra. 
Details of the FUSE instrument have been given bv iMoos et al.l (|2000l ) and iSahnow et all (j2000( ). For the uniformly 
filled MDRS aperture, FU SE achieved a spectral resolution of 0.08 A. Due to thermal misalignments between channels, 
the Jovian spectrum was not recorded in all four channels. Only the LiFl channel, coupled to the Fine Error Sensor 
used for guiding and tracking, shows the expected level of disk emission, decreasing over time as the slit drifted off the 
limb. The disk emission recorded by the SiCl channel has a significantly lower signal/noise, as the effective area of 
this channel is a factor of 4 lower than that of LiFl. The LiF2 and SiC2 channels were off the disk for the entire time, 
showing only emission from the lo plasma torus. The data were reprocessed with the final FUSE pipeline, CalFUSE 
version 3.2. For each day, the data from the first half of each orbit are co-added giving an integration time of 6330 
and 7336 seconds, respectively, for the two days, and the extracted fluxes were converted to average brightness (in 
rayleighs) in the 4" x 20" aperture. 




Figure 1. Snapshot montage of the MDRS slit positions for the FUSE observations of the Jupiter limb, with a time interval of 24 minutes. 



The observations of Saturn were made during 5 contiguous orbits starting at 3:20:06 UT on 2001 January 15. Each 
2000 s exposure was performed using the 30" x 30" slit (LWRS) centered on the planet. At the time of the observation, 
Saturn's heliocentric distance was 9.11 AU, and its geocentric distance was 8.59 AU. The integrated emission from 
the entire Saturn disk was recorded in each observation, with an effective spectral resolution of 0.12 A. The data 
were also reprocessed with CalFUSE version 3.2, and all five orbits of data were co-added together and rebinned to 
0.053 A bins. These observations have also been analyzed bv lGustin et al.l (|2009l ). who showed that the southern polar 
aurora contributed most of the H2 emission through electron impact excitation, using archival HST /STIS FUV-MAMA 
images and spectra taken on 2000 December 7-8. This conclusion is also supported by the narrow line widths, ~0.10 
A, considerably narrower than if the 19."3 Saturnian disk was the primary source of the emission in the 30" x 30" 
slit. The sub-Earth latitude of Saturn was —27.6° so that the southern polar region of Saturn was clearly visible. In 
contrast, both fluorescent and electron impact sources fill the entire aperture in the Jupiter observations. 



2.2. FUSE Spectra of Jupiter and Saturn 
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Figure 2. Selected region from the FUSE spectrum of the 
Jovian limb, with the p hotoclcctron impact excitation model 
HWolven fc FeldmanI 119981 1 overplotted in red. H2 fluorescence 
pumped by solar O VI A1031.91 and Ly/3 produces the excess 
brightness at 1163.8 A, and the feature at 1166.76 A (labeled), 
respectively. 
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Figure 3. Same as Figure [2] for the Saturn disk. The red model 
shows the H2 emission generated by coUisional excitation with au- 
roral electrons. 
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Figure 4. SOHO/SUMER solar Lyman-/3 profiles (from ILemaire et al.|[2002l 'l showing positions of the H2 Lyman (6— O)P(l) line and the 
absorption profile for a column density of lO'^" cm~-^. 



A section of the Jovian spectrum, dominated by electron- impact excited bands of H2 , is sh own in Figure [2l The 
electron impact spectrum was modeled using the method described in lWolven &: FeldmanI ()1998[ ). with a photoelectron 
energy of 30 eV. The spectrum is well described by a model with an H2 temperature of 700 K, and an H2 column of 
1 X 10^*^ cm~^, as shown by the red line in Figure [21 The H2 column density roughly corresponds to the depth in the 
atmosphere where CH4 absorption becomes optically thick. The excess emission of the Q(3) line of the Werner (1—4) 
band at 1163.8 A over the electron excitation model is due to pumping by two solar lines, N III A989.79 and O VI 
A1031.91. The Saturn spectrum, shown in Figure [3l is dominated by the auroral emission, with the corresponding 
electron impact model for Trot = 400 K overplotted in red. 

In both spectra, the most striking feature unaccounted for by electron impact excitation is the P(l) line of the Lyman 
(6—3) band at 1166.8 A, whose position is indicated in the figures. This line is broader, asymmetric, and its peak is 
shifted blueward relative to the electron excited features. The Ly(6— l)P(l) line at 1071.6 A and the Ly(6— 2)P(1) 
line at 1118.6 A have a similar appearance. These lines are part of the fluorescent progression pumped by the solar 
Lyman-/3 line overlapping the Ly(6— O)P(l) line at 1025.93 A. Their shape reflects the fact that there are more solar 
Lyman-/? photons available on the blue side of the Ly(6— O)P(l) line center than on its red side, as can be seen in 
Figure S) This reciprocity in the line profile can only be achieved if we allow the photons in the line wings to scatter 
coherently, preserving their energy offset from the line core. The coherent scattering process requires dropping the 
CRD assumption in the fluorescence model, and using the PRD matrix to compute the emerging line profile. We 
construct a fluorescence model for the H2 molecule, and show that the use of PRD gives a good agreement with the 
data, while the CRD approximation is inconsistent with the observations. Since the electron-impact spectrum has 
been modeled separately, in the following section we will describe in detail the modeling of the fluorescent cascade 
following pumping by solar Lyman-/3. 
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3. FLUORESCENCE MODEL WITH PARTIAL FREQUENCY REDISTRIBUTION 

3.1. Radiative Transfer 

The specific intensity of radiation, 1^, is defined as the energy per unit frequency emitted per unit area in unit 
time, and in unit solid angle along the specified direction, in units of ergs s~^ cm~^ sr"-'^ Hz~^. The emergent intensity 
at the surface of the medium, measured by the observer, is also called surface brightness. The behavior of the specific 
intensity at any point in the medium is described by the radiative transfer equation 

n • V/,y = -K,./^ +>, (1) 

where is the absorption coefficient (cm^^), and jj^ is the cmissivity (ergs s"-'^ cm~^ sr~^ Hz"-'^). This equation 
becomes 

/^-T^ = + j^, (2) 

dz 

for a plane-parallel atmosphere, /i being the cosine of the angle between the direction of the ray and the normal to 
the atmosphere. For every absorber, we can write the absorption coefficient in terms of the scattering cross-section at 
frequency v, as k,^ = nu^, where n is the number density (cm~'^) of the population with cross-section Uv This term 
can represent either line or continuum absorption. In the case of line absorption, n will be replaced by the number 
density of the absorber in the lower state of the transition (n, ). In general, will be the sum of multiple such terms, 
corresponding to all the species and transitions taken into account. Defining the optical depth as 

Ti,{z) = / K^{z')dz' = N(7v, (3) 



Jo 

where N is the column density (cm~^) of the population with cross-section (7,^, and the cross-section is depth- 
independent, the radiative transfer equation becomes 

A^T^ = + (4) 

dTu 

where ~ ju/ Hu is the source function. To simplify the remaining equations, the following derivations are performed 
for /i = 1. The results can be easily generalized for an arbitrary incidence angle. 

The absorption and emission coefficients are the sum of contributions from continuum scattering and individual 
lines. In the case of spectral lines, these coefficients are related to the transition probabilities by 

Ki;j[v) = {hv/4Tr){niBij(f>{v) - njBjiiJj{v)) 
= {hv/4Tr)njA.j^il;{iy) 

where j denotes the upper level, i the lower level, and the angle dependence has been omitted. The quantities nj and 
ni are the number densities of atoms/molecules in the upper and lower level, respectively, and arc related through 
the statistical equilibrium equations. The transition probabilities are given by the Einstein coefficients Aji, while the 
Einstein B coefficients for absorption and stimulated emission are given by 

2 

B^, ^ ^ ^A,,, (6) 

with g^s being the statistical weights of the corresponding levels. 

The radiative transitions are not infinitely narrow at the transition wavelength, but are broadened by the finite 
lifetimes of the excited states, and by bulk molecular motions in the medium. The frequency dependence of the 
absorption and emission line profiles are given by and ipi^)^ respectively. 

The lifetimes of molecular levels, defined as the inverse of Ayiji = i = J2v" ^v'v", determine the natural width of 
the lines, or the Lorentz profile. The gaussian distribution of the molecular motions results in the Doppler profile, with 
the Doppler width given by Ai^jj = vqv/c, where v is the sum of thermal and turbulent motions, v= y/v'^ + b'^. The 

thermal velocity is v^ = y^2kT/m, while b represents the microturbulent velocity, and will be referred as the Doppler 
parameter . The combination of the natural and Doppler broadening results in a profile for the absorption line (j) of 
the form 

H.)^^Hia,u), (7) 

where H(a,u) is the Voigt function, with a = -p-f — and u= ^tt"^. 

In general, the emission profile tp is different from the absorption profile, 0, for the same transition. The two profiles 
arc interrelated through the redistribution function i?(i^,n,j/',n'), that represents the probability distribution that the 
photon is absorbed at frequency v in direction n, and re-emitted at v' in direction n'. It satisfies the normalization 
condition 



(4^) 
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and the emission and absorption profiles can be obtained as 

1 



V'(i'',n') = j R{u,n,v',n')dvdn. 



(9) 



The redistribution functions have been introduced by iHummerl (I1962D for resonance li n es, and have b een gi ven later 
a broa der foundation and extended to subordinate transitions bv iMilkev "eFall (fT975l) : iHeinzell (ITOSl . and IHubenvl 
(|1982l) . The redistribution function describes the coherent scattering of radiation in the rest frame of the emitter (atom 
or molecule), combined with the frequency incoherence determined by the motion of the emitter, the broadening of 
the energy levels, and the collisions with other particles. In this general scenario, the radiation is only partially 
redistributed (PRD) in frequency. The case when this process involves scattering in lines with the same upper level 
but different lower levels (subordinate), is referred to as cross-redistribution (XRD). 

A special case occurs when the absorption and emission frequencies are considered independent, equivalent to the 
complete frequency redistribution (CRD). With this assumption, the redistribution function takes the simple form 



R{v, n, v', n') = 0(i^, n)(j){v' , n') 



(10) 



and the emission profile is identical to the absorption profile. This approximation is the method of choice for small 
optical depths when the Doppler core of the line is predominant, and for forbidden lines occurring between sharp 
energy levels. However, this approximation has important consequences for the photon escape probabilities, and does 
not correctly predict the line profiles. The implications of using PRD vs. CRD in the fluorescence models will be 
further explored in Section fS. 2. 2 1 

Having defined the frequency redistribution functions, the general expression for the source fun ction at frequency 
with multiple overlapping lines and continuum, will be (jUitenbroeklll989l : IHubenv k, Litej|1995() 



S{v) = 



where Qkji{v) is defined as 







k<j 




^^c{i^) + {hv/^'^)Bij(t)i]{v) 


Ui - 


— I + nkQkjt{v) 


Pj J V <Pij{^) 









(11) 



(12) 



with Pj denoting the total probability per unit time for transitions out of level j (radiative and collisional) , and all the 
quantities being angle-integrated. The general form of the angle-averaged redistribution function Rkji{v,v'), taking 
into account collisional broadening, is 



Rkj^iy^v') = (1 - -i)(j)kj{v')(j>tj{v) 



(13) 



where 7 is the coherence fraction, defined as 7 = Pj/{Pj + Qj), Qj being the rate of elastic collisions for level j, and 
Pj is the total depopulation rate of level j. 

We can show that for the H2 transitions pumped by solar Ly/3 in planetary atmospheres, the depopulation rate of level 
j via allowed transitions is much higher than both forbidden and collisional rates. The observed H2 emission is mainly 
coming from the upper parts of the atmosphere, above the altitude at which the methane optical depth becomes one. 
This corresponds to the regions of lowest density and highest temperature in the atmosphere. Computing a weighted 
average of the atmospheric models for Jupiter and Saturn, down to the methane optical depth limit, we obtain H2 
densities of 7x10^^ cm~^ and 10^^ cm~^, and temperatures of 717 K and 350 K for Jupiter and Saturn, respectively. 
Using kinetic theory, the number of collisions per second between H2 molecules is given by Ucoii ~ v^Trd^n-^/SfcT/Trm, 
where d is the diameter of H2 molecules (2.4 A), n is the number density, k is the Boltzmann constant, T is the 
temperature, and m is the mass of the H2 molecules. This expression gives us an estimate of 5x10'^ coUisions/s for 
Jupiter and 5x10^ collisions/s for Saturn, including both elastic and inelastic collisions, between all the hydrogen 
molecules, regardless of their energy or excitation state. This estimate can be considered an upper limit for both the 



elastic collision rate Qj and the inelastic collision rate Cj affecting level j . We can write Pj = Rj 



where i?j is the sum of the radiative rates, Cj the sum of the collisional rates (<10'^ s~^), and Aj the sum of the 
allowed radiative rates (^10^ s~^), and Fj the sum of the forbidden radiative rates (~^10~^ s~^). Since Qj is also 
of the order lO'' s~^, with good approximation Pj ~ Aj ^ Qj. Following this argument, in the fluorescence model 



presented in Section [3.21 we assume 7 = 1. However, for comparison purposes, we also compute the model in the CRD 
approximation, corresponding to 7 = 0. The results for the two limiting cases are shown in Section |4l 



H2 Fluorescence with Partial Frequency Redistribution 
a=0.5 a=1 





a=2 




Figure 5. Frequency redistribution matrices computed for throe values of a, as labeled. The redistribution function is constructed for 
t;D=l km s~^ and a = 0.01, at a resolution of 0.1 Doppler widths. 



In the last term in Eq. [131 Rk]i^^T^') ^he angle-averaged redistribution function derived bv lMilkev et al.l ()1975[) : 
iHubenvl ()1982j ) for transitions with a common upper level and sharp lower levels: 



R''{XI,X2) 



-3/2 



arctan 



arctan 



du. 



(14) 



where xi and X2 are the distances from the line centers of the absorbing and emitting lines, respectively, in units of 



Doppler widths (x = {v—i/q)//S.V£) ). ai is the damping parameter of the absorbing line ( a ; 
and the other quantities are defined as 



y = max{xi — u, a{x2 — u)), u 



y = min(xi + u, a{x2 + u)), a = — 

A2 



\xi - ax2\ 
1 + a 



(15) 



The redistribution matrices depend on the ratio of the central wavelengths of the absorbing and emitting lines, a, 
and on the temperature and turbulent velocities through the Doppler broadening. Therefore, new matrices have to be 
computed for each pair of transitions, and at each depth in the slab, if the Doppler broadening is variable. The number 
of matrix elements to be evaluated is reduced in half by the symmetry property R^^ {xi,X2) = R^^ {—xi, —X2)- The 
form of the redistribution matrices for three values of a is shown in Figure O The color scale stretch in the figure has 
been chosen to emphasize the behavior in the line wings. We distinguish two regimes, one in the line core, where the 
photons are redistributed uniformly forming the Gaussian profile, and another in the line wings, where the photons 
are coherently scattered, preserving the energy difference from the line center. The CRD case can be visualized as the 
central bright square region. It is clear that expanding this behavior beyond the central ±2 Doppler units will result 
in large errors. 

The radiative transfer equation can be solved formally for a medium of total optical depth Tjy, yielding 



Jo 



(16) 



for the reflected and transmitted parts of the intensity, respectively. This expressions can be evaluated directly 
for each frequenc y when the source function a t all the points in the medium is known, usin g either the diff e rentia l 
Feautrier method (jFeautrieiil 1 9641 : [Mihalasll 1 9 78[) . or the integral short characteristics method of lKunasz fc Aueil (jl988[) . 



However, in most cases the source function is itself dependent on the local specific intensity, both through the level 
populations and the emission profiles, as seen from Eq. [TT] 

Even for the cases where the radiative transfer equation can be linearized, the resulting system contains all the 
depth-dependent frequency points and their cross-correlations through the spatial derivatives and the redistribution 
matrices. The direct inversion of such a system is numerically prohibitive. As a consequence, the method of solution 
proceeds iteratively, starting with an initial estimate of the source function to determine the specific intensities and 
level populations, which are in turn used to obtain the new estimate for the source function. This method, known as 
A-iteration, has less memory requirements, but the computing time per iteration is dependent on the number of depth 
and frequency points, and its c onvergence is slow for la rge optical depths. The alternative solution used in this work is 
based on the approximation bv lLiu &: Dalgarnol (I1996D . that docs not require multiple iterations or matrix inversions, 
but relies on a number of simplifying assumptions. This method is presented in detail in the next section. 
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3.2. Fluorescence Model 

The starting model for the treatment of H2 fluoreseenee, described in Section [3.2.1[ was introduced bv lLiu fc Dalgarnol 
(|1996l ). and is based on a set of simphfying assumptions that wiU be used throughout the paper. The model equations 
are derived for a plane-parallel atmosphere, neglecting the stimulated emission term in the absorption coefficient. The 
distribution of molecules in the ground state is assumed to be a Maxwell-Boltzmann distribution given by 

„ (2 J + l)g,exp{~E,j/kT) 

where Ntot is the total column density, E„j is the energy in the {v" , J") level, T is the temperature, and gs is the 
nuclear spin statistical weight. The weight gs is 3 for odd rotational quantum numbers of H2, and 1 otherwise. Unless 
otherwise specified, the quoted temperature is used for deriving both the Doppler line width and the ground state 
distribution. This distribution can be easily replaceable with a user-specified one for special cases. The medium is 
characterized by global temperatures and column densities for all the constituent species. 

The energy difference between the electronic ground state and the excited electronic states is large enough that the 
population of the upper levels is determined solely by the radiative excitation, except when electron impact excitation 
is present. The upper levels of the H2 transitions pumped by solar Ly/3 belong to the excited electronic state B^I]+, 
11 eV above the ground electronic state. At the temperatures in the upper atmosphere of Jupiter and Saturn (300- 
700 K), the average energy of the H2 molecules is 0.03-0.06 eV, insufficient to excite the molecules into the B^E+ state. 
Even though in most parts of the atmosphere the dominant collisions are with molecular hydrogen, these collisions will 
not affect the population of the excited electronic states, because the hydrogen molecules don't have enough energy to 
excite these states, and the collisional frequency between H2 molecules is much smaller than the lifetimes of the excited 
electronic states. Due to the short lifetimes of the B^S]+ state levels ('--^lO^^ seconds), it is unlikely that collisions 
with electrons or hydrogen molecules (of the order of 10^ per second, or less) will change the population of these 
levels before they de-excite into the X^I]+ state via a radiative transition. Excitations and de-excitations between 

neighboring levels in the B^S+ state are dipole- forbidden, and therefore they have a negligible contribution to the 
level populations. The fluorescence rates are also assumed small, so that the distribution in the ground electronic state 
X^E+ satisfies thermal equilibrium, eliminating the need for solving the detailed balance equations. This assumption 
is adequate for planetary atmospheres, but breaks down in the ISM regions where far-UV fluorescence is the dominant 
factor determining the ground state population of H2. Adopting this approximation in general, we can predict the 
effects of PRD on the emergent line shapes, but only approximately the relative line strengths. 

While the photoelectrons (tens of eV, on average) and the auroral electrons (in the 10-100 keV range) have enough 
energy to excite transitions into the B^S+ state, we have taken this contrib ution into account separately, as it's respon- 
sible for most of the far-UV H2 emission observed in planetary atmospheres (jShemanskv fc AiellolllQSStlLiu fc Dalgarnol 
I1996D . The electron-impact excitation model is not restricted to a handful of lines (like the Ly/3 pumping model), but 
includes all measured H2 transitions in the Lyman (B^I]+-X^I]+), Werner (C^n.u-X^I]+), Lyman cascade (EF^I]+- 
BIE+-XIE+), as wen as the B'il]+-Xil]+, D^U^-X^^ , B"^^+-X^j:j, and D'^U.^-X^j:j systems. This electron- 
impact excitation model is shown in Figures [5] and [3] for Jupiter and Saturn, respectively. A similar model for 
continuum-pumped fluorescence, including all far-UV H2 electronic transitions, is especially prohibitive, and we in- 
clude a brief exempliflcation in Section [ 



3.2.1. Single Scattering Approximation 

As can be seen from Eq. I16[ in the absence of emitters in the medium, the emergent intensity through a slab of 
optical depth t^, is I^ir^) = ■ A medium is considered optically thick at frequency v when t(i') ^ 1, and the 

transmitted radiation intensity is attenuated by a factor of e. Let Go{i') denote the speciflc intensity of the incident 
solar radiation field. The total amount of radiation absorbed by the slab will be given by 

g{,^) = goH{i-e~^^''^) (18) 

where t{i>) is the total optical depth summed over all possible transitions and continuum, 

Tii.)^n{iy) + Y,njiiy). (19) 

In this manner, the method accounts for line saturation and line overlaps, and multiple species can contribute to 
T{iy). The frequency dependence of the line optical depth is characterized by the Voigt profile. In the absence of pure 
absorption, photon number conservation requires that the quantity in Eq. I18[ integrated over all frequencies, be the 
same as the total integrated intensity of the fluorescent lines. The fraction of the absorbed intensity from Eq. [18] that 
is due to a particular transition (fc — )■ j) is Tkj{v)/T{v). To obtain the rate at which the upper level j is populated 
through the transition (fc — > j), we multiply this fraction with g{i') and integrate over all frequencies. Denoting 
the emergent intensity at the top of the atmosphere, the total intensity in the fluorescent line (j — > i) will be given by 

^-^d.) ^ (20) 
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where rjj is the prc-dissociation fraction of the upper level, and the last fraction is the branching ratio for the transition 
in line (j — > i) to occur. In molecules, the excited electronic states can overlap with the dissociation continuum of 
another electron ic state, leading to radiationless transitions between these two states, process called pre-dissociation 
(|Herzberglll95nh . This process affects strongly some of the excited electronic states of the hydrogen molecule, such as 
B"iE+, D ^Uu, D'^Tlu, and D" ()Aiello et al.lll984D . A fraction of the population of these states will therefore be 
lost without a radiative de-excitation. This loss is accounted for by the prc-dissociation fraction term in Eq. [20l A 
more in-depth justification of Eg . EOl and its relation to the solution of the radiative transfer equation will be presented 
in Section 13X2] 

This basic treatment of the emission process does not take into account multiple scattering in optically thick lines, 
and cannot reproduce the effects of s elf-absorption. A method to reproduce the self-absorption effects for optically 
thick CO lines has been described in iLupu et al.l (|2007[ ). Due to the integration over frequencies in Eq. [20l all the 
information about the line profile is lost, and the result is equivalent to using the CRD approximation in the source 
function, given that the Doppler core of the lines is not resolved. In the next section we extend this framework to 
allow for depth variations of the model parameters, and to model the emission line profiles using partial frequency 
redistribution. 

3.2.2. Treatment of Frequency Redistribution 

The PRD radiative transfer problem is more complex t han the CRD case, du e to the he avy couplings in frequency 
in addition to spatial correlations and detailed balance (jHubenv fc LitesI 119951 lUitenbroc k 2001). Even restricting 
the geometry to the plane-parallel case, and assuming fixed level populations, the problem remains computationally 
expensive due to intrinsically large optical depths and fine frequency g rids. 

In this section we develop a multilayer extension of the method of iLiu fc Dalgarnol (jl996( ) presented above, that 
reproduces the effects of multiple scattering and the emission line profiles using the frequency redistribution matrix. 
The basic idea is that one uses an equation of the same form as Eq. [201 to calculate the specific intensities at each 
depth in the medium, and then propagate the outgoing intensity to the surface by adding these contribution to the 
next layer, step-by-step. 

We will work in the approximation that stimulated emission is negligible, and that the number density in the upper 
level rij is determined only by the radiative transitions (see also Eq. [SJ. These approximations are justified by the 
short lifetime of the upper level j of the transitions pumped by solar Ly/3. The levels higher than j do not contribute 
significantly to the population of level j via spontaneous emission, since transitions between different ro- vibrational 
levels of the B^I]+ state arc forbidden in the dipole approximation, while the electronic transitions between the B^I]+ 
and X^S+ states are dipole-allowed. The difference between the transition probabilities of the dipole- forbidden and the 
dipole-allowed transitions is about 14 orders of magnitude, and therefore a spontaneous emission out of a ro- vibrational 
level in the B^E+ state is much more likely to go to the ground electronic state X-'^I]+ than to another level within 
the B^E+ state. However, the population of the B^S+ state can also be modified via cascade from the EF^S+ states. 
The EF^I]+ states are populated via electron impact excitation from the ground electronic state, as the transitions 
between these two levels (X and E,F) is also dipole forbidden. This cascade process has been taken into account in the 
electron impact excitation model, computed separately, and added to the fluorescence model for the final spectrum. 

Although the energy difference between the ro- vibrational levels in the B^I]+ state is small enough, such that 
collisional transitions between such levels could be possible given the energy of the collision partners in planetary 
atmospheres, the lifetime of such ro-vibrational levels in the B^I]+ state (~ 10~^ s) is too short for collisional excitation 
and de-excitation to occur. As discussed in Section |3. 11 under the parameters describing planetary atmospheres, the 
upper limit for the average collision rate beween H2 molecules is about 10^"'^ s~^, much smaller than the radiative 
de-excitation rate for a dipole-transition, which is about 10^ s^^. Since this estimate for the collision rate does not 
take into account the kinetic energy of the partners, nor the excitation state of the H2 molecules, the actual collisional 
transition rate between any two excited levels in the B^E+ state will be much smaller. 

Under these approximations, we can write for the population of level j determined by radiative excitation 

k<, J (21) 

k{v) = Kc{v) + {hv/47r)Bij(f)ij{iy)ni. 

After cancellations, the source function will be 

= E£4Sy / Ri^.MH'^ld.', (22) 

^ k<j ' 

where is the general redistribution function (Eq. [T4| for coherent scattering in the atom's frame and negligible 
lower level broadening (metastable state) ()Hubenvlll982l : iHummerll 19621) . A generalization of this expression, including 
the effects of collisions (negligible in our case), is presented in the Appendix. 

The radiative transfer Eq. |4l with the source function given by Eq. I22[ can now be solved using one of the methods 
mentioned in Section [01 We investigated the feasibility of these methods for the full treatment of the II2 molecule. 
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but the large number of transitions (of the order IC*) and the fine frequency and depth grids, renders them impractical 
due to the large mem ory requirements and large number of iterations. On the other hand, the approximation of 
iLiu fc Dalgarnol ()19960 described in Section 13.2.11 is very fast, but suffers from the inability to reproduce the line 
profiles and self-absorption. We are able to show that the PRD process can be easily introduced in this formalism, 
and that the effects of multiple scattering are reproduced by dividing the medium in multiple layers, and applying this 
method layer-by-layer. 

We divide the medium into M layers of constant parameters (number densities, temperatures, etc.), indexed by 
starting with at the surface. The formal solution for the inward radiation (Eq. [T6|) for one layer then becomes 



/-(T'+i,i^) = /-(r',zy)e-^"^'')+e-'^"('^) / e'Sit)dt = r{T\i^)e-^^'''''^ + Id{T^+\i^), (23) 

Jt' 

where Ar = t'+^ — t' and Id(T'+^) represents the local diffuse radiation field. Denoting the incident intensity by 
Goii^'), each layer will be illuminated by go{i'')e~*^'^ \ Substituting this initial value in the the source function from 
Eq. [211 we obtain 



Rl',Miy')e«W)e-'^''>d,''dt, 



or, collecting the optical depth dependent terms. 



-,-At 



-dt{v) 



dv'. 



(24) 



(25) 



Making the usual discretization approximation where all parameters are constant for each layer, will represent 
the number density in level k for the layer I, and k(j/) the absorption coefficient for the same layer. We can also write, 
for short, k{v) = k, n{y') = k' , t{i') = k{i')z = kz, t{v') = k{v')z = n' z. The optical depth integral outlined in the 
square brackets then yields 



K = 



gZ' + ^(K-K') _ gZ'(K-K') 



K=0 



1 



(26) 



Separating the emission process from the absorption process, in the last expression we made the approximation that 
the absorption coefficient will be at the emitted frequency such that the equation for the local intensity becomes 
linear: 



/-(r'+\.) 



k<j 



Ar(i/') 



dv'. 



(27) 



In the last step we introduced the column density in level k for layer Z, ATV^ = UkAz = nfc(z'+^ — z'), and replaced 
the absorption coefficient by the optical depth of the layer (Ar = kAz) in the denominator. 

In the absence of collisions, Pj is the sum of all radiative transitions from level j, containing both the emission 
and absorption rates. In our case, we do not have absorption from the excited electronic states, so that Pj reduces 
to the sum of all emission transition probabilities, or Pj — Aj. Wc also need to multiply this expression with the 
pre-dissociation factor (1 — r/j). In the CRD approx imation, when the redis tribution function takes the form in Eq. llOi 
the expression for intensity Eq. [57] reduces to the ILiu fc Dalgarnol (|1996l ) form in Eq. [201 after taking into account 
Eqs.[n|and[3]for the line optical depth. 

In our multilayer method, the Eq. [571 is only the initial estimate of the intensities at each spatial point. In the next 
step, we start propagating these intensities outward, assuming total reflection at the bottom boundary. In this process 
we account for absorption and re-emission of the outgoing radiation in each layer, writing 



i+i 



-Ar 



"P~^ 4^' 

k<j 



'+1 u') 



Ar(i/') 



.-Ar(i.') 



diy' . 



(28) 



This expression ensures photon number conservation and progressively adds the contribution of each layer to the total 
emergent intensity. The resulting spectrum is computed in a single iteration, and docs not require matrix inversions. 
The computing time depends mostly on the number of grid points in optical depth and the number of transitions for 
which we have to compute the redistribution matrices. 

The properties of this method versus the approximation of ILiu fc Dalgariiol ()1996[ ) (single- layer), and the Feautricr 
solution are shown in Table [T] The output from the three methods is illustrated in Figure [SI for both the PRD and 
CRD cases. This comparison has been performed for a uniform medium of neutral hydrogen, with a total column 
density of lO'^'^ cm^^ and a temperature of 700 K, corresponding to an optical depth at line center of 4x10^. The 
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Table 1 

Characteristics of the radiative transfer methods compared. 





Lambda Iteration 


Multilayer 


Single 




(Feautrier solution) 




layer 


Resources 


High 


Moderate 


Low 


Consistency 


Can be unstable 


Constrained 


Constrained 


Convergence 


Scales writh optical depth 


2-pass layer-by-layer 


Single step 


Spatial variations 


Yes 


Yes 


No 




Figure 6. Comparison of the three methods of calculating the emergent intensity of Ly/3 for an optical depth at line center of 4x10^, 
using the partial redistribution (PRD, left) and the complete redistribution (CRD, right) scenarios. 



frequency grid was centered on Ly/3, with a bin size of 5.44 GHz, and the optical depth grid was spaced using 10 
points per dex at line center. Taking into account the instrumental resolving power, the results from the multilayer 
and Feautrier methods are effectively indistinguishable. The difference in the total integrated brightness of the line is 
about 2% for the PRD case, and less than 10% for CRD, while the computing times differ by a factor of about 100. 
The single layer approximation overpredicts by a large amount the number of photons in the line core, illustrating 
its inability to account for multiple scatterings (self-absorption). The decrease in the Ly/3 intensity seen in the other 
cases is due to its branching ratio of 0.88. 

3.2.3. Effects of PRD Compared to CRD 

As Figure |6] shows, the use of the PRD matrix results in a large effect on the line shape, due to a large number 
of photons escaping in the line wings. This effect becomes significant when the optical depth at line center exceeds 
10^, and is illustrated by the FUSE observations of the atmospheres of Jupiter and Saturn presented in the next 
section. Large variations in the intensity of the exciting radiation over the line profile can also lead to a shift of the 
peak brightness from the line center, when the largest number of photons is absorbed in the line wing and re-emitted 
coherently. 

The use of PRD has even stronger consequences for the predicted integrated line strengths and the peak intensity at 
the line center, as illustrated by Figure [T] The panels show the dependence with total optical depth of the emergent 
number of photons in the central wavelength bin (left) and the integrated number of photons in the line (right), as 
predicted by the PRD and CRD scenarios (in black and red, respectively). The behavior of an optically thick line, 
where the absorption takes place for a set of transitions with a common upper level, is plotted in the upper set of 
panels, while the lower ones display the results for an optically thin line with the same upper level. To guide the 
eye, two vertical lines mark the points where the optical depth is one, and where the column density is lO'^^ cm~^, 
respectively. The transitions are pumped by the solar Ly/3 profile normalized to an integrated intensity of 1. The 
saturation of the fluorescent line brightness at high optical depths in Figure [7] reflects the absorption of all available 
photons from the incident radiation. 

First, we see that the total brightness of the optically thick absorbing line is greatly increased by using PRD, 
due to the fact that more photons are able to escape in the line wings, without being re-emitted in an optically 
thin line from the same progression. In the CRD scenario, more photons are found at the line centre and undergo 
multiple scatterings, increasing the number lost to other transitions. It also shows that in the CRD case the integrated 
intensity remains roughly constant beyond unit optical depth, as expected. Correspondingly, the number of photons 
in the central wavelength bin is greatly suppressed for the optically thin subordinate lines. This fact is determined 
also by the coherent scattering in the line wings, as fewer photons will be redirected towards the line core. The total 
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Figure 7. The dependence of the number of photons in the central wavelength bin (left) and of the total number of photons integrated 
over the line profile (right) on the optical depth for an optically thick line (upper panels) and an optically thin line with the same upper 
level (lower panels). The results from the models using PRD and CRD are shown in black and red, respectively. 



OPTICAL DEPTH =3x10= 

1QQ L 




Figure 8. The dependence of the measured signal (number of photons per resolution element) on the instrumental FWHM. The integration 
using the model with PRD is shown in black, and the one for the CRD case is shown in red. 



intensity is also somewhat decreased in the PRD case, by comparison with CRD, but this effect is not as noticeable 
on individual lines, because the scattered photons are shared among many transitions with the same upper level. The 
peak brightness at the line centre in the optically thick lines remains basically unaffected, reflecting the identity of the 
PRD and CRD matrices in the line core. 

When the PRD treatment is employed, both subordinate (optically thin) and resonant (optically thick) lines show 
lower intensities at the line centre, and the line profiles are widened by the photons escaping in the line wings. This 
process results in a decrease of the line contrast across the spectrum, and is investigated in Section [4. 2[ in relation to 
the search for fluorescent emission from molecular hydrogen in reflection nebulae. Depending on the spectral resolution 
of the instrument, the PRD assumption predicts that the lines will be more difficult to detect than expected based on 
the estimates from CRD calculations. This difference, as shown by Figure [71 will significantly affect the subordinate 
lines. We estimate the effect on the line detection ability by integrating the number of photons from the optically thin 
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Figure 9. The atmospheric models of Jupiter (left) and Saturn (right) used for the synthetic fluorescence spectrum. 

subordinate line in a resolution element equal to the FWHM of the spectral resolution of the instrument. Figure [8] 
shows the resulting number of photons versus resolving power for an optical depth at the wavelength of the absorbing 
line of 3x10'^, for both PRD (black) and CRD (red) models. The model predictions become indistinguishable when 
the instrumental FWHM is larger than lA (R<1000), but we should expect a factor of 5 drop in the signal due to 
PRD when the instrumental FWHM is O.OlA (R=100,000). 

The next section shows specific examples where the use of PRD is required to explain the observations, and the 
multilayer method is used to generate synthetic spectra that reproduce the line profiles for planetary atmospheres. As 
a generalization, we use a toy model to estimate the effects of PRD on the continuum-pumped H2 fluorescence, with 
application to high resolution observations of reflection nebulae. 

4. RESULTS 
4.1. Planetary Atmospheres 

For both planets observed by FUSE we compute the H2 fluorescence spectrum pumped by solar Lyman-/? using 
the multilayer approach with depth dependent parameters. The variation of gas composition and temperature with 
altitude was determined from current atmospheric models. The species with contributing opac ities were H I, H2, an d 
CH4. For the exciting radiation field we used the solar Lyman-/? flux from the compilation of ILemaire et al.l (j2002l) . 
measured on a date close to that of the FUSE observations. 

4.1.1. Jupiter 

The atmospheric model for Jupiter is shown in the left panel of Figure [9l The fluorescence model was integrated 
for a set of 19 lines of sight through the atmosphere, starting with the innermost region falling into the slit, at 13'.'66 
from the planet center, and continuing with a logarithmic spacing up to the outer radius. The depth variation of 
the model parameters was adjusted for each beam, accounting for the variation of layer thicknesses and composition 
with projected distance from the disk center. The observed integrated brightness was also corrected for the slit filling 
fraction during the exposure, as the slit was scanning the limb (see Figure [IJ . The effective area of the disk observed 
is given by A^ff = J A{t)dt, where 

Ait) = 20(4 - <v,) + (i?2^p/2)(asm(10/ii'j„p) - 10/i?,7„p - + sin{d)) (29) 
cos{Q) = {Rju-p - rfo + tYs)/R.jup- (30) 

The Jovian radius Rjup is 20'.'34 (71,492 km), the scanning velocity of the slit v^ is 0'.'00383 s~^, and the initial position 
of the leading edge of the slit, do, is 2'.'68 from the edge of the disk. The integration results in 45% slit coverage during 
the night time exposure. This scaling is employed in extracting the FUSE spectra shown in Figure 1101 Assuming 
that the limb starts approximately 0'.'384 from the outer edge, we obtain that the 83% of the recorded emission was 
originating in the disk, and 17% in the limb. A similar method was used to compute the weights of each pencil beam 
integration to the final brightness predicted by the model, by replacing the planet radius by the position of the pencil 
beam as the limit for slit coverage in each case. This composite PRD model, shown by the red line in Figure [TUl is 
in good agreement with the data, within the error bars. The blue line in the same figure represents the predictions 
of the CRD model, integrated in the same manner, while the green line shows the contribution from the electron 
impact model. Note that photon number conservation is satisfied by both CRD and PRD models. This comparison 
shows that the PRD model is a better description of the data, while the CRD approximation is clearly incompatible. 
Equivalently, the two models correspond to the limiting cases of 7 = 1 and 7 = in Eq. 1131 respectively. In addition 
to the energetic and collision rate arguments, this comparison between the CRD and PRD models and the data is 
further evidence that for the radiative transfer of the H2 lines pumped by solar Ly/3, collisional processes are negligible. 
CoUisional scattering would reduce the amount of coherent re-emission in the line wings and redistribute more photons 
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in the line eore, contrary to observations. The remaining differences between the model and the data could be due to 
unaccounted for variations of the atmospheric composition across the limb of the planet. 



Lyman (6-1) P(1) Lyman (6-2) P(1) Lyman (6-3) P(1) 

, — ' — ' — ' — ' — , — ' — ' — ' , — . — I 1001 — ' — ' — ' — < — ' — ' — ' — ' — — ' — ' — T^n 1001 — ' — • — ' — ' — ' — ' — • — ' — n • 




1071.0 1071.5 1072.0 1118.0 1118.5 1119.0 1166.0 1166.5 . 1167.0 

Wavelength (A) Wavelength (A) Wavelength (A) 



Figure 10. FUSE MDRS Jupiter data, witii PRD (red), CRD (blue), and electron impact (green) models overplottcd. The residuals are 
shown in the bottom pannel, after subtracting both the PRD fluorescence model and the electron impact ecitation model (red and green). 




Figure 11. FUSE LWRS Saturn data, with PRD (red), CRD (blue), and electron impact (green) models overplottcd. The residuals are 
shown in the bottom pannel, after subtracting the sum of the red and green curves. 



4.1.2. Saturn 

As opposed to the the emission generated by electron impact excitation, restricted to the auroral region, the solar 
pumped fluorescence is integrated over the entire disk of the planet. The fluorescence models were calculated for a 
plane-parallel atmosphere, using the depth profil es for the ga s dens ity and temperature shown in the right panel of 
Figure [SI The atmospheric model was derived by iMoses et al.l (j200C)D . A single pencil beam integration was sufficient 
to reproduce the FUSE observations covering the entire disk of the planet. The PRD, CRD, and electron impact 
models are overplotted in red, blue, and green, respectively, on the FUSE spectra in Figure [TTJ The Saturn data have 
a higher signal-to- noise than in the Jupiter case (used the 30"x30" aperture), and we obtain a better agreement with 
the PRD model, while the CRD approximation is again inconsistent with the data. 

4.2. Preliminary Models for Continuum Pumped H2 Fluorescence 

The occurrence of far-UV H2 fluorescence in reflection nebulae can be confirmed by infrared observations of transitions 
from higher ro- vibrational levels [v' > 2) in the ground electronic state, through which the fluorescent molecules cascade 
down to lower states. At expected gas temperatures below 1000 K, the LTE population of these hot ro- vibrational levels 
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would be insufficient to produce the observed infr ared line strengths. Such is the case for the infrared observa tions of 
the reflection nebulae NGC 2023 and NGC 7023 (iGatlev et all 119871: [Martini et al.lll999l: iTakami et al.|[2000l ). which 
support the presence of significant far-UV H2 fluorescence. Observations of these regions have been performed with 
lUE and HUT longward of 1200 A ()Francell2005D . and with FUSE at shorter wavelengths. Fluorescent H2 emission 
has been detected in the lon ger w avelength lUE spectra, namely the feature s at 1575/1608 A, but not conclusively 
in the FUSE data (iFrancd l2005l ). The rocket observations of NGC 2023 (jBurgh et all I2002D also suggest a far- 
UV excess above the dust-scattered light, that might be due to molecular hydrogen fluorescence, but no individual 
transitions have been positively detected. We suggest that for the large optical depths characterizing these objects, 
the coherent scattering of the photons in the line wings is more important than previously thought, leading to a 
decrease of the line contrast. In this case, a high resolution spectrum can look qualitatively different than predicted by 
CRD models. Investigating this effect requires a model of the continuum-pumped H2 fluorescence, including frequency 
cross-redistribution between all excited transitions (> lO'^). Preliminary results show that the expected effects on the 
output spectrum are reproduced by the use of PRD. 
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Figure 12. Comparison between tlie H2 spectrum predicted by the CRD (black) and PRD (red) toy models for NGC 2023. The top 
panel shows a region of the spectrum at the resolution of the numerical computation, while the middle and bottom panels show the same 
models convolved with 0.33 A and 4 A FWHM gaussian kernels, respectively. 

We construc t a toy model, wi th an incident 22,000 K blackbody spectrum normalized to 5000 times the Habing 
galactic mean (jMever et al.ll200ll ). illuminating a uniform II2 cloud with a rotational temperature of 1500 K, a Doppler 
b parameter of 1.8 km s~^, and an integrated column of 5x10^" cm~^. The radiative transfer was performed in a 
single layer approximation, to minimize computational resources. Excitation from higher ro-vibrational levels has not 
been taken into account. The results shown in Figure [12] represent the first PRD calculation of a fiuorescent molecular 
spectrum pumped by a continuum source. The upper panel shows a magnified section of the synthetic spectrum, 
with the PRD model overlaid in red on the CRD model. We note in the PRD model a decrease in the line contrast 
and a blending of the line wings, leading to the emergence of a quasi-continuum, where a significant fraction of the 
scattered photons is found. According to Figure [H we expect that the line peaks would be reduced even more when 
using a multi-layer approximation. The inclusion of dust scattering in the model would have a similar effect, effectively 
increasing the optical depth of the medium. The next panels show the same model convolved with gaussian kernels of 
0.33 A and 4 A FWHM, respectively. These spectra correspond to the resolution of the FUSE and lUE observations, 
and demonstrate that at the resolving power of FUSE there are differences between the two models that could be 
detected, while at the lUE resolution the two models are indistinguishable, and we expect similar predictions for the 
1575/1608 A features. In this sense, this model is qualitatively consistent with both the low signal strength of the 
fluorescent emission in the FUSE data, and the emission features detected by lU E. 
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FU SE observations of planetary atmospheres have revealed the need to include PRD in radiative transfer models 
for H2 fluorescence. As shown in Section [3.2.2[ the inclusion of PRD will affect not only the integrated line strengths, 
and thus the escape probabilities, but also the estimated detectability of the lines due to the spread of the photons 
in the line wings. We constructed synthetic spectra for the atmospheres of Jupiter and Saturn, taking into account 
the depth dependence of the model parameters. The results are in very good agreement with the data, showing 
that the line profiles are reproduced only when accounting for PRD. Using an approximate integration method, we 
show that the use of PRD can significantly affect the resulting spectrum for the continuum pumped H2 fiuorescence, 
containing fewxlO^ overlapping transitions. The result is a decrease in the line contrast that can negatively infiuence 
the detectability of H2 fiuorescence with high resolution spectrographs. 

The H2 lines pumped by solar Ly/3 are very suitable for investigating the effects of PRD on the observed spectrum at 
high spectral resolution, due to the strong variation of the exciting spectrum over the absorbing line profile. Fluorescent 
molecular lines pumped by strong stellar features might also offer the possibility of detecting molecular hydrogen in 
extrasolar planets. In the most favorable circumstances, we find that the maximum ratio of such lines to the stellar 
continuum, assuming that the star and the planet are not spatially resolved, is only a fewx 10~^. Similar considerations 
should be applied to cases where multiple absorbing lines overlap over a narrow wavelength range, but observationally 
they are harder to disentangle. 
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APPENDIX 

This appendix presents the generalization of the source function used in Section 13.2.21 to include the collisional 
contributions. Taking into account that the transition rate out of level j is Pj ~ ^^(C'jfe + Rjk), we can write the 

k 

population of level j as 

".^ = p = p + p = + ' (1) 

^3 ^3 ^3 

where Cij and Rij are the rates of inelastic collisions and radiative transitions, respectively. 

Substituting this expression in the line emission coefficient ( in Eq. [Sj or following the first sum sign in the numerator 
of Eq. [m ), we obtain 

J^3{^) = ^A,,D,,{v)^ (2) 



where we defined 

Dtj{v) = (t^^]{v) 



k 



+ lT,,{u), (3) 



and 



^ ^ / I{^')RU^.^'W. (4) 

k ^ 



On the other hand, the population of level j due to radiative transitions is 

= E / I{^')<t>kA^')dv' = / T,3{^)dv, (5) 

the last equality following from the normalization properties of the redistribution function, Eq. [9l Substituting Eq. [5] 
and Eq. [T] into Eq. |3l we obtain 

D,,{v) = nf<j>,,{iy) + (1 - l)nf<j>,,iiy) + ^T,,{u). (6) 

This expression clearly separates all contribution to the line emission coefficient, and it is easy to implement, given 
that is easy to evaluate from Eq. [5l when the quantities Tij{y) have already been computed. We see that in the 
absence of collisions {Cij = and 7= 1) Dij{v) reduces to Tij{y), and we recover the Eq. [22l after substituting the line 
emission coefficient into the source function. 
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